Unconventional superconductivity in YNi 2 B 2 C 
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We use the semi-classical (Doppler shift) approximation to calculate magnetic field angle- 
dependent density of states and thermal conductivity k zz for a superconductor with a quasi-two- 
dimensional Fermi surface and line nodes along k x = and k y = 0. The results are shown to be in 
good quantitative agreement with experimental results obtained for YM2B2C (Ref. [l|). 

PACS numbers: 74.20.-z, 74.20.Rp, 74.25. Fy, 74.70.Dd 



I. INTRODUCTION 

YM2B2C is a type II superconductor with a relatively 
high transition temperature T c — 15.6K. 2 Although ini- 
tially thought to be a conventional s-wave supercon- 
ductor, accumulated evidence soon suggested otherwise. 
Power law behaviour in the heat capacity C p /T oc T 2 was 
the first indication that YN12B2C is an unconventional 
superconductor with point nodes in the gap function!^ 
However, the field dependence of the heat capacity was 
found to be C p oc \/H, indicative of line nodes4>£ The 
NMR spin relaxation rate \jT\ was measured to be oc T 3 
with no Hebel-Slichter peak£, again consistent with line 
nodes in the gap function. Finally, Raman scattering 
showed a peak in the electronic A\ g and i?2g response/ 
possibly indicative of a Bi g symmetry gap function. Such 
a gap function takes the form A(fc) oc k x k y , and thus 
has symmetry required line nodes along k x = and 
k y = 0.— ^ In contrast to these findings, field-angle depen- 
dent measurements of thermal conductivity and specific 
heat were claimed to be indicative of point nodes in the 
gap function^ ' 10 ! 11 The reconciliation of these results, 
and hence the symmetry of the gap function, remains an 
important unresolved issue. 

YM2B2C belongs to the crystallographic space group 
I4/mmm (No. 139, D^ h )H The lattice is body-centred 
tetragonal (bet) with a = 3.5261 and c = 10.543A 3 
According to symmetry analysis for crystals*^ gap 
functions with line nodes are found only for singlet pair- 
ing, while point nodes are found only for triplet pairing. 
Various nodal configurations can occur, depending on the 
irreducible representation of D^h by which the supercon- 
ducting order parameter transforms. Nodes in the gap 
function are normally detected via quasiparticles (q.p.'s) 
which appear in the vicinity of gap nodes in fc-space as a 
result of either finite temperature, impurities, or Doppler 
shift in the presence of an applied magnetic field. In these 
kinds of measurements, the nodes will be invisible if there 
is no Fermi surface in the direction of the nodes, thus the 
shape and connectivity of the Fermi surface plays an im- 
portant role. 

The Fermi level crosses the 17th, 18th and 19th bands. 
The topology of the Fermi surface is highly sensitive 
to the precise position of the Fermi level due to a dis- 
persionless band between the T and X points. Thus 



different band structure calculations share common fea- 
tures but the resulting Fermi surfaces have significant 
differences! 13 ' 14 i 15 Yamaguchi et al^ correlated their re- 
sults with de Haas-van Alphen (dHvA) measurements in 
order to fix the Fermi energy. The 18th and 19th bands 
were found to produce closed Fermi surfaces around var- 
ious points in the Brillioun zone, however the 17th band 
produces a large electron Fermi surface multiply con- 
nected by necks. Part of this surface appears as dHvA 
oscillations perpendicular to the c-axis. The orbits do 
not appear to be closed in the c direction; instead they 
seem to possess a two-dimensional character that extends 
in the c direction, as evidenced by the upward curvature 
of the dHvA frequencies about the [001] direction, shown 
in Fig. 3 of Ref. IS 

In the vortex phase of a type II superconductor q.p.'s 
may be either localised about vortex cores or delocalised. 
It was shown some time ago that the contribution to 
the low-energy density of states in a superconductor with 
line nodes comes from delocalised q.p. 's in the vicinity of 
the nodes,— The delocalised q.p.'s can be treated with 
a semi-classical (Doppler shift) approximation; this ap- 
proach provided a good description of field-dependent 
specific heat and thermal conductivity of the line node 
superconductors YBa 2 Cu30 [17( and CeCoIn 5 Q. How- 
ever, Volovik's argument does not extend to point node 
superconductors. For point node superconductors, the 
semi-classical calculation may still be performed^ but, 
as may be expected, these results are not in agreement 
with any experiment involving putative point node su- 
perconductors so far. 

In this article, we use the semi-classical approximation 
to calculate the field-angle dependent density of states 
and thermal conductivity for a superconductor with line 
nodes and a quasi-2D Fermi surface, for the purpose of 
demonstrating that the results of such measurements on 
YM2B2C are in fact consistent with this scenario, in con- 
trast to what has been claimed i 

For simplicity, we assume that the Fermi surface has 
the shape shown in Fig. [TJ for which the q.p. energy 
spectrum takes the form 
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FIG. 1: (Colour online) Fermi surfaces for three different val- 
ues of b = — . The surface in c) has the largest value of 

b. 



metry has line nodes along k x = and k y = 0. The 
Fermi momenta along these nodes (parametrised by k z ) 
are 

&fi,3 = (0, ±(2m(e_p — e' F cosckz)) 1 ^ 2 , k' F sincfc 2 ) 

w (0,±fe F ,fe^sincfe z ) (2) 

fcF2,4 = (±(2m(ei? — coscfcz)) 1 / 2 , 0, &^sincfc z ) 

~ (±kp, 0, k' F sin ck z ). (3) 

The magnetic field rotates in the xy plane with an angle 
e with respect to the x axis, 



H = iJ(cose,sine,0). 



(4) 



The supercurrent circulates perpendicular to the field as 
a function of the distance r from the vortex core and 
winding angle /?, 

v s( r ) = r — (— sine cos j3, cose cos j3, sin /3). (5) 
2mr 

The Doppler shifts associated with each line node are 
ati(r) = k Fl ■ v s , 
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[±k F cos e cos /3 + &^ sin /3 sin ck z \6) 



ot2,4( r ) = = n [t^f sin e cos (3 + k' F sin /? sin cfc z ](7) 



II. DENSITY OF STATES 

In the semi-classical treatment, the argument of the 
Green's function iu n is replaced by ioj n + a where a is 
the Doppler shift. The quasiparticle energy is E(k) — 



^Je 2 {k) + A 2 (fc) which is « y vf^fc 2 + w 2 fc| in the vicin- 
ity of a nodeJ^ Here fci points in the direction of the 



node, /c2 is perpendicular to k% in the xy plane, and the 
gap velocity is v g — SA/c^lnode- In the vicinity of the 
jth node, the Green's function takes the form 



G(k,iu n ,r) 



iuj n + OLj{r) + vpk\ 
{iuo n + ctj(r)) 2 + vpk\ + v 2 g k\ 



(8) 



where iuj n = iuj n + iTq and To is the scattering rate at 
zero energy. The density of states is 



N(u, r) = V 3G(fc, Q, r). 

7T * J 



(9) 



We divide the volume of integration into four curved 
cylinder-shaped volumes, each centred around a line node 
on the Fermi surface 19 and perform the integration across 
the disk spanned by k\ and hi 
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where po is the integration cut-off. In the clean limit 
|aj/r | ^> 1 the density of states is 



N(0,r) 
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(11) 



Averaging over the vortex cross-section, we obtain 
1 1 



(N(0,r)) H 



drr / dp dk z 
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4ir 3 v F v g irR 2 Je Q 
x(|oi(r)| + \a 2 {r)\ + \a 3 (r)\ + \a 4 (r)\). 

(12) 



This leads to the result 
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where b = v F /v F , c is the lattice constant in the c direc- 
tion, C = cose, S = sine and E is the complete elliptic 
integral of the second kind. Using N F ~ l/cv g £o and 
£o/R ~ JHJTT7 2 , we find 
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FIG. 2: Clean limit density of states (the dimensionless ex- 
pression in the square brackets of Eq. 1 14 p as a function of 
rotating field angle e (in radians) for b = 0.02, 0.05 and 0.10. 



where ks is the Boltzmann constant and v F is the Fermi 
velocity in the direction of k. By again dividing the vol- 
ume of integration into fo ur regions and introducing the 
integration variable p = » /vpk 2 + v 2 k 2 we find 
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Using ((HE]), in zero magnetic field we get 
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This function is shown in Fig.[2]for various values of b. It 
is seen that deviations from the perfectly 2D cylindrical 
Fermi surface leads to a softening of the cusps in the 
density of states. 

In the dirty limit jct^/Tol < 1 we get 
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In a finite magnetic field, terms linear in the Dopplcr 
shift will vanish upon integration. So the magnetic part 
of the thermal conductivity is 
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which produces no oscillations with respect to the rotat- 
ing field. 
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III. THERMAL CONDUCTIVITY 



In the clean limit, this reduces to 



The thermal conductivity tensor is given by the Kubo 
formula, which is expressed in terms of the imaginary 
part of the Green's function as 
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(16) The integrand is 
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The off-diagonal components vanish in the vortex aver- Again, the off-diagonal components vanish in the vortex 
age, and the diagonal components are average, and the diagonal elements are 
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(25) 



where K is the complete elliptic integral of the first kind. 
k zz is plotted in Fig. [3] for different values of b. In the 
limit b — > the cusps are sharp, however the oscilla- 
tion amplitude goes to zero. The oscillation amplitude 
increases exponentially with b. 
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FIG. 3: Clean limit thermal conductivity n zz (the dimensional 
expression in the round brackets of Eq. I25|l as a function of 
rotating field angle e (in radians) for b = 0.02, 0.05 and 0.10. 
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In the dirty limit, ([20]) reduces to 
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Similar to the dirty limit density of states (fTS")) . there are 
no rotating field-dependent oscillations in the dirty limit 

of K zz . 



IV. DISCUSSION AND CONCLUSIONS 



The topolo gy o f the true Fermi surface of YM2B2C 
shown in Ref. [l5j is difficult to discern, however the va- 
lidity of our calculation only requires that the Fermi sur- 
face exists at the positions of the nodes and spans all 
or most of the Brillioun zone in the c direction with a 
slight curvature characterised by the parameter b. The 
main point is that the cusp features observed in the 
field angle-dependent heat capacity-^ are a feature of 
line nodes and the cusp features observed in the field 
angle-dependent thermal conductivity are a feature of 
line nodes on a quasi-2D Fermi surface. In contrast, the 
semi-classical (Doppler shift) treatment of point nodes 
(which may not even be valid) is insensitive to Fermi 
surface topology and produces neither cusp features nor 
four-fold oscillations 



Using R « 4 x 1CT 8 m (for a 1 T field), T 

15 „,„ -j q w in7 



0.56 



K, Ep — 9 Ryr^ v F w 3 x 10 7 m/s, gap maximum 



1 K leads to 

3 ~ ]^Q4 



(26) 



An = Vghkp = 30 K and scattering rate Tq 
an estimate of the prefactor in (|2"5|) of <^ -^B- Ie. 
W/Km. The experimentally observed oscillation ampli- 
tude is » 2 x 10 -3 . Comparing with the oscillation ampli- 
tudes shown in Fig. [31 one may deduce that the value of 
b is approximately 0.02. Such a small value of b produces 
sharp cusps in the field angle-dependent k zz oscillations 
and is therefore fully consistent with experiment. 

Thus the most straight-forward model that best de- 
scribes accumulated observations on YM2B2C is that the 
superconducting order parameter is A ~ k x k y , which be- 
longs to the irreducible representation Bi g of the point 
group D^h-, with associated line nodes along k x = and 
k v = 0. 
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